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We present a theoretical analysis of the selective darkening method for implementing quantum 
controlled-NOT (CNOT) gates. This method, which we recently proposed and demonstrated, con- 
sists of driving two transversely-coupled quantum bits (qubits) with a driving field that is resonant 
with one of the two qubits. For specific relative amplitudes and phases of the driving field felt by the 
two qubits, one of the two transitions in the degenerate pair is darkened, or in other words, becomes 
forbidden by effective selection rules. At these driving conditions, the evolution of the two-qubit 
state realizes a CNOT gate. The gate speed is found to be limited only by the coupling energy J, 
which is the fundamental speed limit for any entangling gate. Numerical simulations show that at 
gate speeds corresponding to 0.48J and 0.07J, the gate fidelity is 99% and 99.99%, respectively, 
and increases further for lower gate speeds. In addition, the effect of higher-lying energy levels 
and weak anharmonicity is studied, as well as the scalability of the method to systems of multiple 
qubits. We conclude that in all these respects this method is competitive with existing schemes for 
creating entanglement, with the added advantages of being applicable for qubits operating at fixed 
frequencies (either by design or for exploitation of coherence sweet-spots) and having the simplicity 
of microwave-only operation. 
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I. INTRODUCTION 



An important step towards the realisation of a quantum computer is to implement a set of universal gates from 
which all other operations can be composed. The quantum controlled-NOT (CNOT) gate is a prominent example of 
a two-qubit universal gate, requiring only the addition of single-qubit gates to form a complete universal set. The 
physical implementation of a gate depends on the type of qubit that is used, but even for a specific qubit type a 
certain gate can often be implemented using a variety of methods. The ideal method delivers a gate that is fast 
compared to qubit decoherence times, has a high fidelity, is simple to implement and does not introduce constraints 
that compromise the coherence time of the qubits or the implementation of other gates. 

One classification of two-qubit gates is based on the required nature of the interaction between the qubits pQ. 
When the interaction term contains diagonal elements in the single-qubit energy eigenbasis, often referred to as 
longitudinal or zz-coupling, the transition frequency of one qubit depends on the state of the other qubit. Although 
this spectroscopic shift enables simple resonant driving for all operations |2J 13] , the shift also leads to continuously 
evolving phases that must be compensated by refocusing schemes In contrast, for transverse coupling the energy 
splitting of one qubit does not depend on the state of the other qubit. In this respect, the system can be described 
as a set of effectively uncoupled qubits [FJ |S] . The desired entangling evolution can then be induced by dynamic 
manipulation of the system. Clearly, the degeneracy due to the lack of a spectroscopic splitting forbids any method 
based on frequency-selectivity. Focusing on superconducting qubits [7Hl5j. previous methods used either additional 
coupling elements [T6H29] , extra qubit states outside the computational basis (25J I30H32] or shifting levels in and out 
of resonance by DC [33H35] or strong AC fields [3M3H]- 

Another important consideration is that a method can require certain bias conditions of the qubits, which might 
not always be compatible with optimal coherence time conditions, or the implementations of other gates. This is 
especially important for qubits with a so-called sweet-spot, an optimal bias point at which the qubit is rendered 
insensitive to specific noise channels, and coherence times can improve by orders of magnitude. For some qubit types 
and coupling schemes these sweet-spots are intrinsically connected to having a transverse coupling. Previous work 
[5] l6[ |2"T1 l36[ I39MI] has specifically focused on two-qubit quantum gates that can be fully implemented at sweet-spots. 

In this work we theoretically analyse a CNOT gate based on the selective darkening method as proposed and 
demonstrated in reference [39 . This method is for transversely-coupled qubits and requires driving two qubits 
simultaneously with a single frequency, with specific relative amplitudes and phases. The method was developed in 
the context of superconducting qubits, but can also be used for any other coupled-qubit system with pure transverse 
coupling [32] • The basic principle of the gate is explained in section [nj for both small and large coupling energies. 
In section [III] we briefly discuss the relation of selective darkening to other comparable methods. The time-domain 



evolution of the gate is investigated in section IV including the possible gate errors caused by ac-Stark shifts due to 
off-resonant driving of other transitions in the system. A numerical study of the gate evolution, gate errors and gate 
speed is presented in section [V] The influence of possible higher- lying levels is analysed for weakly anharmonic qubits 
in section |VI| Lastly, the scalability of the method to systems of multiple qubits is studied in section |VII[ and we 
finish with the conclusions and discussion of the results. 



II. SELECTIVE DARKENING 



A. Basic principle 



Selective darkening provides an entangling operation between two coupled qubits. To operate the gate, one drives 
both qubits with a common frequency, but individually tuned amplitudes and phases. In this section we first give 
an intuitive explanation of the method, and for simplicity assume a coupling energy between the qubits that is much 
smaller than the other energies in the system. 

The selective darkening method is suitable for any coupled-qubit system where the effective coupling term is purely 
transverse, i.e. where the matrix describing the interaction, written in the energy eigenbasis of the uncoupled system, 
does not contain any diagonal elements. Here we consider a system of two transversely coupled qubits that is described 
by the Hamiltonian 

H a = -\{^+^)+Ja^af\ (1) 

where A,; is the single-qubit energy splitting of qubit i, J is the qubit-qubit coupling energy and o x %,z are the 
Pauli spin matrices (cri^.z = ^x,y,z <S> 1, &x)y,z = 1 ® &x,y,z)- Note however that other (effective) Hamiltonian with 
pure transverse coupling, i.e. with any combination of ai^yd'^l coupling terms, would not change the principle of the 



3 



method. The same remark holds for systems where the coupling is not direct between the two qubits, but mediated by 
an extra coupling element such as a harmonic oscillator [43] or a SQUID 22 (Superconducting QUantum Interference 
Device). In figures 1(a) and 1(b) two relevant example systems are depicted. 




FIG. 1. Selective darkening (a) Schematic diagram of two inductively coupled flux qubits and two antennas for individual 
driving, (b) Schematic diagram of two transmon or charge qubits, coupled via a harmonic oscillator, and with two antennas for 
individual driving, (c) Energy level diagram of the coupled-qubit system. Arrows of the same color indicate transitions of the 
same qubit and are degenerate in frequency, (d) The normalized transition strengths of the four transitions in (c) as a function 
of the ratio of driving amplitudes ai/(ai + 02) for ip2 — <pi = 0. For </?2 — pi = n the dashed and solid lines are interchanged. 
The black dotted lines indicate the condition for which the corresponding transition is darkened. 



As a first step we treat the coupling strength J as a perturbation, taking J <C |Ai — A2I. We also assume that 
I Ax — A2I <C Ax, A 2 , so that we can ignore terms containing the ratio J/Aj. With no loss of generality, we also take 
Ax > A2. To first order in perturbation theory, the four energy eigenstates are given by: 

|0) = |00) 

|1> = |01) - — — ^— — |10> 
Ax — A 2 

|2)H10) + A7^|01) 

|3) = |11). (2) 

The energies are unaffected by the perturbation to first order. In other words E% — Eq = E3 — E\ = Ax and 
Ei — Eq = E3 — E 2 = A 2 . (In the next section it is shown that the relations E 2 — E a = E3 — E\ and Ei—Eq = E3 — E 2 
hold to any order). The energy levels are shown schematically in figure 1(c). The arrows indicate the transitions of 
interest; the blue and red arrows describe the transitions of qubit 1 and 2, respectively. As mentioned above, each pair 
of transitions is degenerate, which is typical for transverse coupling |44j . Let us now add a driving term of frequency 
ui, resonant with one of the degenerate pairs of transitions, that couples with amplitude ai and phase (pi to qubit 1, 
and a 2 and ip 2 to qubit 2: 

-Hdrivc = ai cos(wi + (fii)^^ + a 2 cos(cot + (p 2 )a^ 

= ^^ + 6^-^, (3) 



where 



^ ± ri vc = ye ± ^ai 1 ) + |e ± -ai 2 ). (4) 



The amplitudes ai and a 2 are real and positive, and the phases ipi and ip 2 are real. If we choose u) = A 2 /h (in 
our first-order approximation), we match the transition frequency for flipping the state of the second qubit, i.e. we 
should drive the transitions |0) f-> |1) and |2) «-» |3) [corresponding to red arrows in figure 1(c)]. We now evaluate the 
transition strengths — (l \ -f/drive |&) of both transitions, i.e. the matrix elements that govern the Rabi frequencies 
of the oscillations (ujr = 2\T\/h). The term -ffdrivc represents the (time-independent) co-rotating field, meaning that 
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we take the rotating wave approximation, giving -ffdrivc = H^ rivc for fc > Z and Hdrivc = H drivc for k < I (The 
derivation of the co- and counter-rotating terms is done in appendix [A]). The resulting transition strengths are: 

<1| ffdrive |0> = -YA^Az^ 1 + (5a) 

(3| H dTive |2) = +^^-L_ e -^ + |e"^. (5b) 

Crucially, the two transition strengths are not equal. Figure 1(d) shows the normalized \T\ — \T\/(a\ + 02) as a 
function of a\j (a\ + 02), for a fixed phase difference tf2 — <Pi = 0. The blue lines correspond to the two transitions of 
qubit 2. The difference in transition strength leads to a different evolution for the two corresponding transitions. This 
can be conveniently visualized using a Bloch sphere representation, as shown in figure 2(a). The black and grey arrows 
both represent the state of the qubit that is resonantly driven (in our example qubit 2), where black (grey) indicates 
that the other qubit is in state |0) (|1)). The transition strength is a complex number and both the amplitude and 
phase are relevant for the evolution of the state. Difference in the amplitudes of the two transition strengths leads to 
different Rabi frequencies of the oscillations (u/r^i ^ uir^)- Difference in the phases of the transition strengths leads 
to different rotation axes (</?_r,i 7^ fR,2)- 

We make some simple observations. A single-qubit gate is a gate that changes the state of one qubit with the 
change being independent of the state of the other qubit (s). In the Bloch sphere picture this means that when a 
rotation is induced on one qubit, the rotation speed and the orientation of the rotation axis do not depend on the 
other qubit being in state |0) or |1). This situation is indicated in figure 2(d). This only occurs when the driving 
field couples exclusively to the qubit that is resonant with the driving field; in our case this corresponds to a\ = (or 
&i/(ai +02) = 0, see figure 1(c)). All other settings of the driving field lead to entangling evolution of the system, and 
can be used to create entangling gates. Note that if one intends to perform a single-qubit gate, and the drive is not 
applied exclusively to that qubit, the two-qubit evolution that is aimed for in this work should be taken into account 
as a possible source of errors, unless the system has a tunable coupling that can be reduced to zero [TBI l24l 145] . 



(a) general (b) selective darkening (c) cross resonance (d) single-qubit gate 

|o> |o> ' |o> |o> 




ID ID ID ID 



FIG. 2. Bloch sphere representation of the evolution of a resonantly driven degenerate transition for different driving conditions. 
The black and grey arrows show the evolution of the state of the target qubit, when the control qubit is in state |0) and state 
|1), respectively, (a) General case. The two state vectors rotate with two different Rabi frequencies, and around two different 
axes, (b) Selective darkening. The rotation of one of the state vectors is fully suppressed. The other state- vector has a non-zero 
Rabi frequency, (c) Cross resonance. The state vectors rotate around the same axis, with equal Rabi frequency, but in opposite 
directions, (d) Single-qubit gate. The state vectors rotate around the same axis, with equal Rabi frequency and in the same 
direction. 



Even though many different settings of the driving field can be used to create entangling gates, two special cases 
stand out. The first is where the driving field resonant with one qubit is exclusively applied to the non-resonant qubit 
(a 2 = 0, ai/(ai + a 2 ) = 1). This driving scenario will lead to equal rotation speed of both state vectors, but opposite 
direction. This case will be discussed further in section HTT1 

The second special case is the main focus of this paper, where one of the transition strengths equals zero [depicted in 
figure 1(d) by the black dotted lines]. Here one transition is fully suppressed, i.e. darkened, while the other transition 
has a finite transition strength. Specifically, if we choose: 

fl2 J ra\ 
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and (fi2 — <pi — 0, we find that the Rabi oscillation frequency for the transition |0) ■<-» |1) is zero, so that we only 
drive the transition |2) O |3). This condition provides exactly what is required for a CNOT gate: for the appropriate 
driving duration the state of the second qubit is flipped if the first qubit is in state |1), while it remains unaffected if 
the first qubit is in state |0). The 0-CNOT gate, where the second qubit is flipped only if the first qubit is in state 
|0), is achieved by taking equation ^ and ipi — <ys 2 = ft- 

The speed of the gate is determined by the transition strength of the non-darkened transition, which is now given 
by: 

(3|g d rivc|2) =a 2 =ai . 3 . ■ (7) 
Ai — A 2 

One must take care that both a\ and a 2 are small enough to not excite other transitions by off-resonant excitation. 
As a simple estimation for the maximum gate speed we take the driving amplitudes such that for all the off-resonant 
transitions the transition strength is smaller than the detuning of the transition frequency with the driving field: 
(l\ -ffdrivc \k) < \fiio — huji^k\- In this case the transitions |0) O |2) and |1) ■<-» |3) provide the tightest restrictions 
((2| iJdrivc |0) s» (3| i/drivc |1) ~ ai/2), giving ai/2 < (Ai — A 2 ). The resulting estimate for the maximum Rabi 
frequency of the CNOT transition is 

WR, max ~ iJ/h. (8) 

This means that the maximum gate speed is determined by J, which is the fundamental upper limit on performing a 
two-qubit gate; a two-qubit gate cannot be faster than the coupling strength in the system [46]. How close the gate 
can get to this estimate for the maximum speed (for a given fidelity) is studied numerically in section |v} 



B. Strong coupling 

In the simplified calculations of the previous subsection, the energy eigenstates are easily identifiable as being almost 
equal to the states of the computational basis. However, the Hamiltonian [equation ([I])] is simple enough that it can 
be diagonalized without approximations: 

|0) = cos(0i/2) 1 00) - sin((V2) |H) 
|1) = cos(0 2 /2) |01) - sin(0 2 /2) |10) 
|2) = cos(6» 2 /2) |10) + sin(<9 2 /2) |01) 

|3> = cos(6» 1 /2) |11) + sin(0i/2) |00) , (9) 

where 

2J 2J 

tan01 = A7TA^' tan02 = AT^- (10) 



For compactness of the following equations, we also define 

0i + O2 



We now find that 



and 



The condition for selective darkening of qubit 2 transitions is then given by 

a 2 sin 9 + 



(11) 



(l| J ffdrive|O)=-ye-^ 1 sm0 + + ^e-^cos0_ (12a) 

(3| i/drive |2) = +Y e ^ Vl Sin6 + + COS0 - (12b) 



(2| i/drivc |0) = ^ e"^ 1 cos 9+ + %^ sin 6_ (12c) 
(3| ifdovc |1) - ye-^ 1 cos0 + - ^-e-^ sm6_. (12d) 



a 1 



(13) 



G 



with ip 2 - ifi = for the 1-controlled CNOT gate (1-CNOT), and ip 2 - <fi = tt for the 0-CNOT gate. This condition 
generalizes equation ([6]), which was derived for the weak-coupling limit. The condition for selective darkening of qubit 
1 transitions is 

02 = cos 6+ 
ai sin 6- ' 

with 93 2 — fi = 7r for the 1-CNOT gate, and ^2 — <Pi — for the 0-CNOT gate. In this case (i.e., with arbitrary values 
of Ai, A2 and J), the transition frequencies are given by: 

E 1 - E = E 3 - E 2 = l - (V4J 2 + (Ax + A 2 )2 - + (A, - A 2 ) 2 ) , (15a) 

E 2 - E = E 3 - E 1 = i ( + (A x + A 2 ) 2 + ^4J 2 + (A x - A 2 )a) . (15b) 

Note that this derivation shows that the selective darkening gate can also be used when Ax = A 2 , giving a! = a 2 for 
the required amplitude ratio for both qubit 1 or qubit 2 as the target qubit. In this special case the selective darkening 
method is equivalent to the method proposed by Beige et al. |47j and has recently been demonstrated experimentally 

ESI. 



III. RELATION TO OTHER GATES 



Independently of our work on selective darkening (SD), a cross resonance (CR) scheme for realizing two-qubit gates 
for transversely-coupled qubits was proposed in reference |40j . The CR coupling scheme relies on driving one qubit 
(the control qubit) at the resonance frequency of the other qubit (the target qubit) and has the result of driving 
Rabi oscillations in the target qubit in one of two opposite directions based on the state of the control qubit. Even 
though at first sight the driving condition and resulting dynamics might seem to be very different from those of the 
SD scheme, the two schemes are related by a simple rotation. This relation can be seen by looking at figures 2(b) 
and 2(c) and noting that in both cases the two-qubit dynamics involve rotations of the target-qubit state about a 
single axis. As a result, the two-qubit evolution commutes with any single-qubit rotation about the same axis. One 
can therefore convert any CR-like operation into an SD-like operation by applying a single-qubit rotation designed 
to cancel one of the two possible CR rotations. Crucially, because the single-qubit and two-qubit rotations commute, 
they can be applied simultaneously. Starting from the CR scheme as a reference point, one could say that the driving 
field applied to the target qubit in the SD scheme is designed to oppose one of the two rotations in the CR scheme, 
thus producing CNOT-gate dynamics. 

It is interesting in this context to look further back into the history of coupling schemes for superconducting qubits. 
After the proposal of the FLICFORQ coupling scheme [5], it was shown that that this scheme is a special case of a 
family of coupling schemes using the physics of double resonance [SB] [37]. In particular, even though the coupling 
scheme used in reference [55J might seem very different from FLICFORQ, they are both special cases of the double- 
resonance family of coupling schemes. The general condition that needs to be satisfied for double resonance is given 
by: 



± yj {Tiloi - Ai) 2 +aj±yj {huj 2 - A 2 ) 2 + a\ = hui - huj 2 . (16) 

It is then interesting to note that the condition for the CR HO] and SD gates [55], i.e. driving both qubits at 
the frequency of one of them with small amplitudes, satisfies the double-resonance condition. One can then wonder 
whether these coupling schemes can also be seen as two more (closely-related) members of the double-resonance family. 

The CR and SD schemes can indeed be seen as special cases of double resonance. However, they also involve some 
qualitative differences with previously studied cases, such as the FLICFORQ and ac-Stark-shift-induced resonance, 
making them more like distant relatives of the conventional double-resonance schemes (It is worth noting here that 
another different scheme of the same family was discovered in numerical results in reference |37j). One difference 
between the more conventional double-resonance schemes and the CR and SD schemes becomes apparent if one plots 
the energy-level diagrams of the two driven qubits in the dressed-state picture. In the schemes of references [H[55], one 
transition between dressed states of qubit 1 becomes resonant with one transition between dressed states of qubit 2. 
As a result, a swap-like operation takes place with the two single-qubit transitions driven in opposite directions, such 
that the total energy is conserved. In the case of CR and SD, two pairs of transitions are resonant and are involved 
in the gate operation, as can be seen in figures 1 and 4 of reference 0D]. The result turns out to be a conditional 
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operation: the well-known CNOT-gate dynamics in the SD case, and a CNOT gate combined with a single-qubit 
rotation in the CR case. 

The classification of different gates can also be done using the geometric representation of references [JJ HPJ [SO] . 
Being related by simple single-qubit gates, the CR and SD CNOT gates correspond to exactly the same dynamics of 
the Makhlin parameters. The more conventional double-resonance class of gates, as well as the gate of reference |51j . 
all exhibit oscillations corresponding to the iSWAP gate. As a result the Makhlin parameters exhibit similar behaviour 
for these cases. 

In reference [S2] a scheme is discussed to control the individual states in degenerate subspaces of three qubits 
coupled in a loop. That scheme is based on the same principles as the SD and CR method. The gate proposed in 
reference [ST] is similar to the CR and SD method in the sense that it also provides a microwave-only entangling gate 
between two qubits and does not involve shifting dressed states into resonance by strong driving. This gate however 
requires a coupling type that is at least partially longitudinal. Creating entangling evolution between two coupled 
qubits by tuning the amplitudes and phases of a driving field is also discussed in references [53j and |54j . 

Another intimately related effect to SD is illustrated by the dark states encountered in circuit quantum electrody- 
namics. When two qubits are both coupled to a harmonic oscillator, and the qubits are biased to have equal splitting, 
the energy levels exhibit an anti-crossing. Some of the transitions are typically found to be forbidden. This effect can 
be interpreted as SD, as the condition for selective darkening for two qubits of equal splitting is when both qubits 



are driven with equal amplitude [58] (see section II B I . Recently, it has been proposed and demonstrated that these 
properties can be exploited further by defining a new logical qubit from two conventional qubits, using only two out 
of the four original states [55 -57 . Within this definition, the coupling to the resonator can be treated as tunable, and 



the qubit itself can be shielded from decay through the resonator. 



IV. TIME-DOMAIN EVOLUTION OF THE GATE AND SPURIOUS DRIVING OF OTHER 

TRANSITIONS 



To compare different methods for creating entangling gates, or to design the optimum parameters for a certain 
experiment, it is important to know the potential sources of errors, and how these errors depend on the various 
parameters in the system. The first potential source of errors that we study is related to the ac-Stark effect, i.e. the 
energy shift that the driving field can induce on the non-resonant transitions in the system. Specifically, the applied 
driving field is resonant with the qubit that is to be flipped (the target qubit of the CNOT gate) , but the same field 
also drives the control qubit, which typically has a different transition frequency. This off-resonant driving leads to a 
predictable, but unwanted deviation from the pure CNOT-type evolution of the system. 



A. Ideal gate evolution under the rotating wave approximation 



The time evolution of the system is best studied in a suitable rotating frame. In this case the driving term -ffdrive m 
equation |3| can be conveniently split into a (time-independent) co-rotating term V and a (time-dependent) counter- 



rotating term W (see appendix LAI). For the co-rotating term V we find the matrix elements given in equation (12 i-d) 

with -f/drivc = V. 

We now make the rotating wave approximation: The counter-rotating term W oscillates very fast (at frequency 
2ui) and is therefore approximated by its average value, zero. The gate we would like to perform is the usual CNOT 
gate, with qubit 2 as the target qubit. We therefore choose 

hoj = E$ — E2 = Ei — Eq 

= l - (V(Ai + A 2 )2+4/2 _ V /(A 1 -A 2 )2+4J2) (17a) 

* = S ^± (17b) 

<2l COS0- 

¥2-<Pi= 0. (17c) 



We are now left with 



H = ^(Ai-A 2 )2 + 4J2 



-1 











-1 











1 












(18) 



(1|V|0) = 
(2|F|0) = ^ 



(cos 9+ + sin 9 + tan 6-)e 



(3|F|2) =aisin0+e 
ai 



(3|F|1) 



(cos 0+ — sin 9 + tan ) e 



(i\V\j) = (j\V\i) 



(19) 



In a first approximation, we ignore the matrix element V20 and V31 (and their Hermitian conjugates). We then find 
the effective Hamiltonian 



/0 \ 


V> 



:S2 



(20) 



\o v 32 / 



These two terms commute, and their effects can be easily superimposed. If the Hamiltonian is allowed to operate for 
a time t = /j./ (4| Vj3 2 |) = h/{Aa\ sin6> + ), the first term produces the transformation 



/ e V(Ai-A 2 )2+4J2t/2 

e V(Ai-A 2 )2+4J=t/2 


\ 



e ^i v /(Ai-A 2 ) 2 +4J 2 i/2 

e -V(Ai-A 2 ) 2 +4J 2 t/2 



(21) 



and the second term produces the transformation 



(1 
1 





\0 -ie lipi 



-le 



(22) 



These transformations are given in the rotating frame. They can be transformed back to the lab frame to give: 

/ e~ iEot \ 

e- lElt 

e~ lE2t 

\ e~ iE3t J 

The leftmost matrix represents the Larmor precession of the system in the lab frame. Taking <pi = 0, the rightmost 
matrix is exactly the CNOT gate, up to a single-qubit phase gate and a global phase factor (meaning that the Makhlin 
parameters [49] are equal to those of the standard CNOT gate). 



/ 1 

1 



\0 -ie ltpi 



(23) 



B. Single- and two-qubit errors caused by the ac-Stark effect 



We now include the two (or four including the conjugate terms) matrix elements that we have ignored in the 
previous subsection: 

(2| V |0) = y (cos6»+ + sin 6+ tan6L) e^ 1 (24a) 

(3\V |1) = y (cos (9+ -sinfl+tanMe - '* 1 - ( 24b ) 
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These matrix elements produce an ac-Stark shift on qubit 1 (the control qubit). If 9 + and 9- are small, the above two 
matrix elements are approximately equal to a\/2 (up to a phase factor). This matrix element then causes a shift o f 
a\/2(E 2 —E a — Tiuj). This frequency shift can be incorp orate d easily into the transformation given in equations ( A3|A7 |, 



simply modifying the energies that enter in equation ( Al I 



Note however that the two matrix elements above are slightly different, so that the ac-Stark shift of qubit 1 
depends on the state of qubit 2. This situation would induce an unwanted controllcd-phase (CPHASE) gate. If we 
take A 1 /h = 6 GHz, A 2 /h = 5 GHz and J/h = 0.1 GHz, we find 6 1 = 0.0058tt and 2 = 0.063tt, which gives the ratio 
between the matrix elements: 

(g|VjO) = cos g + + sin g + tang, = 1 Q2 
(31 V II) cos 6 + - sin 9+ tan 6L 

This means that for the given parameters, the two different values for the ac-Stark shift will differ by about 4%. 
This 4% difference of the frequency shift should therefore be designed such that it causes an overall phase shift much 
smaller than 2-k over the gate time of typically a few nanoseconds. It can then be neglected. Otherwise, it could 
harm the fidelity of the gate. Even for extremely strong driving (a\/h =1 GHz), which gives a large ac-Stark shift 
(about 0.5 GHz), over a gate time of 4 ns, we obtain an unwanted CPHASE rotation by an angle of 0.167T. For weaker 
driving this unwanted angle will be smaller. 

Note that in our case, this CPHASE-type error can also be corrected relatively easily using single-qubit gates. In 
this case the SD-CNOT part of the gate is made intentionally to deviate slightly from a perfect CNOT gate by using 
a longer or shorter gate time, such that the evolution of the state-vector overshoots or undershoots the normal 7r- 
rotation. If the correction to the gate time is chosen appropriately, the combined evolution of the unwanted CPHASE 
gate and the intentionally driven CNOT gate can be transformed into a perfect CNOT gate using single-qubit gates 
before and after the two-qubit driving. This effect is related to the well-known fact that a CPHASE gate can be 
transformed into a CNOT gate using only single-qubit gates |58| . This type of correction will be demonstrated in 
section |Vj 

A more rigorous analysis (see appendix [B]) , where the driving field is treated quantum mechanically, confirms the 
single-qubit and two-qubit errors found here. In addition, even though we mostly look at the high driving limit where 
the classical description of the field gives accurate results, some phenomena can be identified more easily. Specifically 
it is shown that for increasing amplitudes of the driving field not only the levels of the control qubit are shifted, but 
also the transitions of the target qubit are no longer perfectly degenerate. In this case no transition can be fully 
darkened. The result is a single-qubit-type error on the target qubit. As mentioned above, single-qubit errors can 
be corrected relatively easily. It should be noted however that this error can slow down the CNOT gate. A perfect 
SD-CNOT gate requires the state vectors in the example of figure [2] to be on opposite sides of the Bloch sphere at the 
end of the gate, and co-rotation of the transition that should be darkened leads to a saturation of the effective gate 
speed. This will also be considered in section |V} 



V. NUMERICAL SIMULATION OF GATE DYNAMICS FOR TWO-LEVEL QUBITS 

To quantify the fidelity of the selective darkening CNOT (SD-CNOT) gate, the evolution of the system was simulated 
numerically. We take the driving conditions derived from the simple weak driving case given in section!!!} the frequency 



of the driving field was chosen to be uj = (E± — Eq) /h = (E3 — E 2 ) /Ti, and the amplitude ratio as given in equation ( 13 1. 
The amplitudes are ramped up and down using the half-sine envelope function sin(7rf/i ga te), with t running from 
to the total gate time £ ga te- Note that we chose to gradually turn on and off the driving field rather than using step 
functions, because in many cases the steep edges can lead to undesired excitations in the system. This approach 
increases the gate time with a factor of 7r/2 compared to the minimum gate time where the driving amplitude is 
turned on and off abruptly. Other envelope shapes would give similar correction factors. 

For a given maximum amplitude a\ 1 the evolution of the system was simulated. The evolution operator at the end 
of the gate was evaluated according to its overlap with the ideal CNOT gate, which is quantified using the fidelity 
T = Tr(Mg ato McNOT)/4. Here Af gatc is the (unitary) matrix of the simulated gate and Mcnot the ideal CNOT 
gate from equation p3| ) (which includes the phase factors i). This calculation is performed for 20 different gate times 
around an estimated optimum. An interpolating polynomial function was used to fit the dependence of the fidelity on 
the gate time, providing the final accurate estimate for the optimum gate time. The gate was recalculated using this 
gate time, giving the final result for the fidelity. This sequence is repeated for different driving amplitudes a\ (with 
02 adjusted to the ratio given in equation (|l~3"|). The used qubit splittings are A\/h = 6 GHz and A 2 /h = 5 GHz 
and the phases of the driving fields are (pi = ip 2 — 0. 

The results of the simulations are shown in figure [3] Panels (a) and (b) were calculated for a coupling strength 
J/h = 100 MHz, and panels (c) and (d) for J/h = 25 MHz. The blue circles give the result for the procedure as 
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FIG. 3. Numerical simulation of a SD-CNOT gate for 2-level qubits. Panels (a) and (b) are calculated for coupling strengths 
of J/h — 100 MHz and (c),(d) for J/h — 25 MHz. Blue circles represent the results for the SD-CNOT gate, incorporating only 
corrections for Larmor precession and ac-Stark shift of the control qubit. The other two data sets represent the SD-CNOT 
gate including error corrections using single-qubit phase gates (black diamonds) and arbitrary single-qubit gates (red squares) . 
In panels (a) and (c) the gate error 1 — T is plotted versus the driving strength. Panels (b) and (d) give the corresponding 
gate speed, expressed as the effective coupling strength J e g divided by the natural coupling strength, for which the fidelity T 
is achieved. The black lines indicate the expected dependence from equation (121. The values used for the qubit splittings are 
Ai/ft = 6 GHz and A 2 /h = 5 GHz. 



described above, using the lowest number of corrections. We only take into account the Larmor precession, and for the 
control qubit the ac-Stark and Bloch-Siegert shift due to the (off-resonant) driving, disregarding for the level shifts any 
effects related to the coupling. The bare SD-CNOT gate is calculated first, and the corrections are applied afterward 
as phase gates. The required phase for the ac-Stark and Bloch-Siegert shift corrections is calculated analytically from 
the energy shifts £-E ac _stark = a\/2{E 2 -E - hu) and ^^Bloch-Siegert = a\/2(E 2 -E + huj) and the gate time £ gate , 
taking into account the envelope shape. The correction for the Larmor precession is given by the conjugate of the 
leftmost matrix in equation ( 23 1 . Note however that all of these corrections often do not have to be implemented 
physically in an experiment. The internal phase evolution of a microwave-source that is set to the resonance frequency 
will automatically follow the Larmor precession of the qubit, and a single-qubit phase gate can often be incorporated 
easily in preceding or following single-qubit gates. 

For the black diamonds and red squares we performed the same procedure, but in addition we allowed single-qubit 
phase gates (black diamonds) or arbitrary single-qubit rotations (red squares) before and after the two-qubit gate. For 
each simulation of the two-qubit gate an optimization algorithm is used to find single-qubit gates that optimize the 
fidelity of the gate. This approach allows for correcting certain error types, such as the CPHASE-type error discussed 
in section IV B Note that in all the three cases, the corrections influence not only the maximum fidelity of the gate, 
but also for which i ga t e the maximum is found. 

We characterized the gate speed as the effective coupling strength J cff = h/t gatc ■ n/2, where the factor tt/2 takes 
into account the envelope function of the driving field. In panels (b) and (d), the black lines represent the expected 
gate speed as given in equation (12). For low driving strengths the gate speed follows exactly the linear behaviour 
as expected from the transition strength. The gate error (given by 1 — J 7 ), increases monotonically with driving 
strength, except for the case with arbitrary single-qubit corrections where the error is ~ 10~ 5 throughout the whole 
range and is limited only by the accuracy of the numerical calculation. The low error in the latter case shows that 
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the errors in the SD-CNOT gate are either single-qubit errors, or two-qubit errors that can be incorporated in the 
combined evolution to form a perfect CNOT gate, such as the CPHASE-type errors discussed in section IV B Also 
in experiment this would be an attractive way to enhance the gate performance. 

The slow-down of the gate for the results with arbitrary single-qubit gates (red squares) at higher driving frequencies 
can be explained with the error types discussed in the previous sections. If the transition |0) «-> 1 1) is not fully darkened, 
and the state vector rotates in the same direction as the |2) -n- |3) (CNOT) transition, a longer gate time is required 
to achieve the maximum separation of the two state vectors (to achieve a perfect CNOT gate the alignment of the 
two state vectors should be exactly opposite). Alternatively, a CPHASE-type error can require a slight overshoot of 
the CNOT gate evolution to achieve the best fidelity when including the single-qubit gates. Both processes are likely 
to contribute to the observed slow-down of the gate. 

At higher driving strengths the gate speed starts to deviate from the linear dependence of the weak-driving limit. 
Especially where the error is larger than 0.1, the gate fidelity becomes less well-behaved. For J = 25 MHz, the 
effective coupling strength even seems to exceed the natural coupling strength, but this can be explained simply by 
the fact that the performed gate is no longer close to the desired gate. Also for J = 25 MHz, we see an improvement 
when using the single-qubit phase corrections (black diamonds). 

Even without the extra gate optimizations (blue circles), we find a high fidelity at relatively high gate speeds. For 
the coupling strength of J = 100 MHz, gate speeds corresponding to 0.48 J and 0.07J give a gate fidelity of 99% and 
99.99% respectively. The exact numbers will depend on the system parameters. Note that these results were obtained 
using a fixed driving frequency and amplitude ratio, obtained from simple analytic equations, and using a fixed pulse 
shape. Using more advanced optimization schemes, for example using the DRAG (Derivative Removal by Adiabetic 
Gate) method (59) [60], these number can likely be enhanced further. 



VI. WEAKLY ANHARMONIC QUBITS 



The above analysis of the SD-CNOT gate and the possible errors is accurate for true two-level qubits. However, 
many qubit implementations are in fact multi- level systems. In this case a nonlinear level structure allows the system 
to be used as an effective two-level system, where only two out of a manifold of levels are used as the (computational) 
qubit states. When dealing with weakly anharmonic qubits, i.e., qubits where the transition frequencies between 
consecutive energy levels differ only slightly, the additional energy levels have to be taken into account, as they can 
significantly influence the dynamics of the system. For example, one obvious complication that arises in this case is 
the leakage to these additional states. 

We take the computational basis to be formed by the lowest two energy states of the qubits, labeled |00), |01), 
1 10) and 1 1 1) , and the states outside the computational basis are labeled |02), |20), etc. Taking into account the 
third level of each qubit, and making the approximation of weak coupling (J <C Ai — A 2 ) and weak anharmonicity 
(Sj = hcoi^2:j — ^o-H-iy Ai — A 2 ; note that under this definition 5j is negative for phase and transmon qubits and 
positive for capacitively shunted flux qubits), we find the following approximate expressions for the energy eigenstates: 
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The respective approximate energies are given by: 
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(27) 



Note that the labels are ordered for increasing energy of the states, except for the interchange of labels 3 and 4. This 
interchange makes that states 0-3 correspond to the two-level qubit case, and makes it easier to compare the two 
cases. 

Let us consider the CNOT gate with qubit 2 being the target qubit. In the qubit basis, this means that the |0) — » |1) 
transition is to be darkened, and the |2) — > |3) transition provides the CNOT gate dynamics. The required relation 
between the driving amplitudes is again given by equation ([6]) and (p 2 ~ fi = 0. For simplicity in the expressions we 
also take </?i = y 2 = 0. The transition strengths of the darkened and CNOT-gate transition are now given by 
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(29) 



From these equations we can recognise the first consequence of the additional levels: The transition strength of the 
|3) |2) (CNOT) transition is reduced by a potentially significant factor. This reduction is a consequence of the 

mixing between the states |11) and |20), which introduces a new term into the matrix element (3| ^a^ + ^foa, 2 '' |2). 
For systems where both the anharmonicity and the coupling energy are small compared to the detuning between 
the qubits, the speed of the gate will be greatly suppressed. So in order to achieve large gate speeds with weakly 
anharmonic qubits, one must design the system with one or both of these two factors relatively large. 

Next we analyze the leakage to the states outside the computational basis. The largest leakage arises from the 
transitions |1) — > |4) and |3) — > |6), because qubit 1 (the control qubit) is driven the strongest. The corresponding 
matrix elements and energy detunings are given by: 
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and 



(E 4 - Ex) - (E, - E ) *5 2 - {Ai _l 2)2 (31a) 

372 

(E 6 -E 3 )- (E 3 -E 2 ) = S 2 - — . (31b) 

Ai — A 2 

To avoid leakage, the transition strengths of the unwan ted t ransitions should be lower (possibly much lower) than the 
detuning with the driving fields (similar as in section 



have approximately equal energy splitting (« S 2 )■ The 
Using this transition we find the estimate for the maximum allowed Rabi frequency of the CNOT transition: 



II A I : (^| i/drivc \k) < \huj — hu!i^k\- Both leakage transitions 
3) -H- 1 6) transition clearly has the highest transition strength. 
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(32) 



How close one can get to this estimated maximum is simulated numerically, and will be discussed shortly. 

A third issue that one needs to keep in mind when dealing with weakly anharmonic qubits is that the additional 
energy levels cause frequency shifts that break the symmetry E$ — E 2 = E\ — Eq. In particular, the energies of the 
state 4 (« |02)) and 5 (« 1 20) ) will be close to the energy of the state 3 (« |H))- Fortunately, they cause shifts in 
opposite directions and almost cancel. The net shift is given by (see the expressions for the energies above): 

SF = 2J 2 2J 2 ^ 2,P{8 1+ 5 2 ) 

3 A x — A 2 - S 2 Ai - A 2 + <5i ~ (Ai - A 2 ) 2 ' 1 ' 

This frequency shift is not specific for the selective darkening method, but typical for weakly anharmonic qubits, and 
results in an always-on interaction that drives a CPHASE gate. 

We again perform numerical simulations to determine the total fidelity of all the combined errors. First we start 
with simulations where we vary the anharmonicity in order to see its effect on the gate speed and fidelity. We 
use three different combinations of anharmonicity: (8\/h — —400 MHz, 5 2 /h — —360 MHz), (8\/h — —200 MHz, 
S 2 /h = -180 MHz) and (5 x /h = -100 MHz, S 2 /h = -90 MHz). The other qubit parameters are A 1 /h = 6 GHz, 
A 2 /h = 5 GHz and J/h = 25 MHz. 

The results are shown in figure [4j At low driving strength, the gate speeds have a good correspondence with the 



expected dependence from the transition strength [equation ( 28b |], confirming the slowdown of the SD-CNOT gate for 
decreasing anharmonicity. This implies that for weakly anharmonic qubits it is recommended to design the detuning 
between the qubits to be relatively small, so the negative impact on the gate speed due to the weak anharmonicity 
can be minimized. Due to the low gate speeds, all the energy shifts induced by the driving and the higher-lying levels 
are also more significant. For this reason, the simple analytical calculation of the AC-Stark and Bloch-Siegert shifts 
of the control qubit (as was used for the calculation of the data represented by the blue circles in figure [3]) no longer 
gives useful results. Therefore we do not plot the corresponding results in figure [4] 

Another prominent feature is the decrease of the gate speed at high driving amplitudes. One possible explanation 
for this effect is a competition between the SD-CNOT gate and an AC-Stark induced CPHASE gate. For figure (4^d) 
a plateau is observed at low driving strength. This plateau is related to the CPHASE evolution caused by level shifts 
induced by the higher energy levels of the qubits, and is therefore independent of the driving amplitude. As mentioned 
before, a CPHASE-type error can be corrected by altering the duration of the two-qubit driving and using single-qubit 
rotations. The plateau indicates that the CPHASE contribution to the final CNOT gate is even dominating over the 
selective-darkening contribution. The origin of the plateaus in (d) and (f) at high driving power is not known. The 
same holds for the sharp sudden increase in the error. Note however that at these points the driving strength is 
already higher than would be advisable given the anharmonicity of the qubits. 

A common approach to improve gate speed and fidelity for weakly anharmonic qubits is to use pulse shaping (see 
for example the DRAG method [5H1 [SO]) • This approach could lead to improved gate speed and fidelity, but will not 
be discussed here. 

In the experiment described in [3T], a cross-resonance (CR) gate is performed. Given the close relationship between 
the SD and CR method, we calculate, for comparison, the SD performance for the described system. The qubit 
parameters in this case are A 1 /h = 5.854 GHz, A 2 /h = 5.528 GHz, Si/h = 225 MHz, S 2 /h = 255 MHz and 
J/h = 6 MHz. The coupling J is estimated from the numbers provided in the reference. Figure [5] shows the result of 
the simulations. The gate speed in the simulations is in good agreement with the typical gate speeds in the experiment. 
With these system parameters one achieves a relatively high effective coupling, primarily thanks to the small detuning 
between the qubits. At very low driving strength we see a similar plateau as in the previous parameter set, and again 
attribute this to level shifts due to the states outside the computational basis. From the fidelity data in panel (a), we 
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FIG. 4. Numerical simulation of a SD-CNOT gate for 3-level qubits with weak anharmonicity. Black diamonds represent 
the SD-CNOT gate including single-qubit phase gate corrections and the red squares represent the SD-CNOT gate including 
arbitrary single-qubit gates. Panels (a,b) are calculated for anharmonicities of 8i/h = —400 MHz, 62 /h = —360 MHz, (c,d) 
for 5i/h = -200 MHz, 5 2 /h = -180 MHz and (e,f) for Si/h = -100 MHz, 8 2 /h = -90 MHz. In panels (a,c,e) the gate error 
1 — T is plotted versus the driving strength a\ . Panels (b,d,f ) give the co rresp onding gate speed J e s / J for which the fidelity is 
achieved. Dashed lines indicate the expected dependence from equation ( |28b[ |. The qubit parameters used are Ai/h = 6 GHz, 
A2//1 = 5 GHz and J/h = 25 MHz. 



conclude that the leakage to these states is negligible at these driving strengths: if there were significant leakage to 
these levels, the results using the arbitrary single-qubit corrections represented by the red squares could not lead to 
the gate improvement that is observed, because the single-qubit corrections that we apply in the calculations cannot 
transfer back any state-population that is lost during the SD gate. The increase of the gate error for the results with 
arbitrary single-qubit gate corrections at ai/h = 0.1 GHz is not yet understood. Also we do not see the significant 
slow-down of the gate at high driving strength as was measured in the experiment. A possible explanation for this is 
that different pulse-shapes were used, that in the experiment also DRAG-type 59, 60 pulses were applied, or that the 
presence of the resonator has an additional influence that is not covered by our direct-coupling model. Nevertheless, 
the calculation illustrates that similar performance can be expected for the SD and CR methods. 



VII. SCALABILITY 



Let us now consider the scalability of the SD-CNOT gate to systems of more than two qubits. We take a system 
composed of a one-dimensional chain of N transversely-coupled qubits and use the Jordan- Wigner (JW) transforma- 
tion to turn it into a system of fermions occupying a one-dimensional lattice with N sites. The transformation is 
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FIG. 5. Numerical simulation of a SD-CNOT gate for 3-level qubits with parameters similar to the experiment in reference [41] . 
The qubit parameters used are Ai/h = 5.854 GHz, A 2 /h = 5.528 GHz, 5i/h = 225 MHz, 8 2 /h = 255 MHz and J/h = 6 MHz. 
Black diamonds represent the SD-CNOT gate including single-qubit phase gate corrections and the red squares represent the 
SD-CNOT gate including arbitrary single-qubit gates. In panel (a) the gate error 1 — J 7 is plotted versus the driving strength 
ai . Panel (b) gives the corresponding gate speed J B ft / J for which the fidelity is achieved. The black line indicates the expected 
dependence from equation (28b I. 



given by: 
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where c and are annihilation and creator operators, 

/ n n \ 

(35) 

and \ + _ = l' 1 ' ® l' 2 ' • • ■ ® 6% g + _ ■ • ■ ® 1^. This transformation transforms our original Hamiltonian, i.e. 

H = E -Y*® + E (36) 

into 

AT iV-1 

# = E ^i^tci + E J M+i (4^+1 + cJ +1 Ci + cjcj+i + Ci+iCj) + const. (37) 

»=1 i=l 

This Hamiltonian does not look much simpler than the original one. In order to obtain something simple out of it, 
we now write the Hamiltonian as 

* / + + \ / i ,\ transpose 

H = [c\, ■ ■ ■ ,ff N ,c 1 ■ ■ ■ ,c N J M [ci, ■ ■ ■ ,c N ,c\, ■ ■ ■ ,& N ) + const., (38) 

where M is a 27V x 27V matrix. It is always possible to define new operators 7, (with i — 1,2,..., TV) that also obey 
Fermi statistics and for which the Hamiltonian takes the simple form 

„/, , \ „ / . . \ transpose 

H = {%, ■ ■ ■ ,7^,71, • • • ,7jvJ M (71,- • • ,7at,7i, • • • ,T N J + const., (39) 

with a diagonal matrix M (note that this step can be carried out by diagonalizing the matrix M). In other words, 

N 

H = Y^ Eajji + const. (40) 
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The physical meaning of the above Hamiltonian is the following. We have defined new creation and annihilation 
operators, which are combinations of the creation and annihilation operators for the original fermions. With the new 
definition, the Hamiltonian has turned into a Hamiltonian of non- interacting fermions. By diagonalizing the matrix 
M, we therefore obtain the available single-particle energy levels for these newly defined particles. This transformation 
can be interpreted as going from a description using the bare qubit states to a picture where the states are dressed 
by the qubit-qubit coupling energies. For a system with fixed qubit-qubit couplings, the dressed states are the most 
natural and physically relevant states to describe the dynamics |61j . 

The total energy of the entire system is obtained by identifying which single-particle energy levels are occupied and 
taking the sum of the energies of these levels. In other words, any energy eigenstate can be expressed as 7] 7] ■ ■ ■ T k \G), 
where k denote the occupied single-particle states and \G) is the ground state of the system. This derivation 

provides a somewhat intuitive explanation to the result of reference |62j , that the Larmor frequency of a given qubit 
is independent of the states of the other qubits. Since the original Hamiltonian is equivalent to a Hamiltonian that 
describes non-interacting fermions, it is natural that the change in energy induced by the addition of one particle 
into a certain energy level is independent of whether the other energy levels are occupied or not. This statement 
is independent of the strength of the coupling: even for strong coupling, where the energy eigenstates of the spin 
system are very involved (which corresponds to the fermion energy eigenstates having a large extension in space), 
the symmetric energy-level structure is still preserved. The only requirement is a one-dimensional chain with a 



Hamiltonian of the form given in equation ( 36 1 



Now let us turn to the eigenstates of the Hamiltonian. In general, determining these eigenstates is harder than 
determining the eigenvalues, because even the ground state has a complicated form. The ground state \G) is defined 
by the set of conditions 

7, \G) = (41) 

for all values of i. In other words, there are no 7 fermions in the ground state. However, in the language of the 
c fermions, the above set of equations are not transparent at all and do not have any simple solution in general. 
This difficulty could make calculations complicated. However, we are mostly interested in the practical situation 
where the typical scale of the coupling energies (J) is small compared to the typical energy scale of the difference 
between different qubit splittings {5 A). We can therefore perform perturbation-theory calculations where we start by 
considering the case J = first and then analyse the effect of adding a small coupling term. 

In the case of zero coupling (J = 0), the matrix M is already diagonal, which means that the 7 fermions coincide 
with the c fermions. This means that each (single-particle) energy level is localized at one lattice site with no tunneling 
between these sites. The state of the entire system is defined simply by specifying the occupation (0 or 1) of the 
different lattice sites. We now add a small coupling term (J 7^ 0). For simplicity we first consider the condition that 
the coupling strength is very small compared to the qubit Larmor frequencies (J <C A; note that here we use A, not 



5 A) and make the rotating wave approximation, meaning that the terms in equation (37) that create or annihilate 
two c fermions are ignored: 



N N-l 



H « ^2 AiC^i + ^2 Ji > i+1 + + const - ( 42 ) 

»=1 i=l 

This approximation leads to a conceptual simplification, because we now know that creating a 7 fermion in a given 
energy level corresponds to creating a c fermion in a quantum state that is (in general) spread over the entire 
lattice, with no involved mixing between creation and annihilation operators. Another important result of the above 
approximation is that the ground state reduces to a very simple form. Regardless of the exact form of the operators 
ji in terms of the operators £j, the set of conditions 7$ \G) =0 leads to the set of conditions Cj \G) = 0, such that the 
ground state corresponds to zero occupation of all the single-particle states in both languages. Note that the number 
of particles is also the same for any state in both languages &c — 7^7) • 

With the above approximation, we have turned the Hamiltonian into a Hamiltonian that conserves the number of 
particles (even in the language of the c operators) . Note also that there are no interaction terms in the Hamiltonian. 
We can therefore proceed in calculating the form of the operators 7$ by considering a single-particle problem (This is 
allowed because in the absence of interactions the form of the different single-particle wave functions will be insensitive 
to whether other wave functions are occupied or not). For this purpose we imagine that there is only one particle in 
the system, and this particle can hop (i.e. tunnel) between the different lattice sites. 

When calculating the form of the single-particle wave functions, we use the condition J -C SA. This condition 
implies that although the c fermions can hop between the different lattice sites, the hopping matrix elements are 
small compared to the detuning between neighbouring sites. As a result, the different single-particle states (i.e., those 
described by the 7 fermions) will each be concentrated around one lattice site (recall that when J = the states are 
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completely localized at individual lattice sites). One can calculate the spread of the states using perturbation theory. 
Writing 



N 



7> 

we find that: 

i/>i(i) = 1 + 



J 

SA 



2 



/ /• i o\ J»,*±lJi±l,i±2 > n( J 

^ (l±2)= (A i -A i±1 )(A i -A i±2 ) +0 U 
/ j \ I* — J I 

■0*0') ~ ( J when |i - j| » 1 

; (44) 

Note that the above relation can be inverted straightforwardly: 

N 



Ci 



(45) 



In principle the localization situation changes when two or more lattice sites are degenerate. For example, in the 
case of two-fold degeneracy and an otherwise symmetric environment, the (single-particle) energy eigenstates will 
be symmetric and antisymmetric superpositions of the particle being at one of the two sites. However, provided 
that these degeneracies occur only for largely separated lattice sites, the effective tunneling matrix element will be 
exponentially small in the distance between the two sites in question: (assuming j > i) 

^effective Aj — A i+1 Aj - A J+2 Aj - Aj_i ^ 

We can therefore ignore such long-range tunneling and focus on the practically relevant case where the 7 fermions 
and the c fermions are almost equal, with exponentially decaying tails describing the spread of the 7 fermions. This 
situation is related to the phenomenon of Anderson localization, where disorder causes the single-particle states to be 
localized with exponentially decaying tails. 

We can now start doing the calculations for the Rabi frequencies. Before going into long expressions, we make the 
following observation. If a driving signal is applied to qubit i, the driving term in the Hamiltonian can be expressed 
in the language of fermions as: 

tfdrivc = ]J i 1 - 2£ ] £ j) ( 5 * + 4) COs(L0t). (47) 

j<i 

The relevant matrix elements for purposes of evaluating the Rabi frequency of qubit k therefore have the form: 



(0| 7m In ■■■% 



j<i 



^•••TiClO), (48) 



where the index i represents the driven qubit and the operators on the far right and far left (i.e. those with indices 
m,n,...,p) describe the state of the surrounding qubits. Let us now assume that i and m are separated by a large 
distance. From the localization argument above, we know that the fermion created by the operator 7^ is concentrated 
mostly on one side of site i (i.e. to the right if m > i and to the left if m < i), with only a very small probability of 
occupying site i or crossing to the opposite side. As mentioned above, this probability decreases exponentially with 
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the distance between i and m. We can therefore make the approximation that the operators (ci + c|J and 7^ anti- 

commutc and that the operators Ylj<i fl — 2 4^) an< ^ Tm either commute or anti-commute, depending on whether 

m > i or m < i. Note here that the operator Yij<i fl — ^ctcyj coun t s the number of particles to the left of site i and 
produces a sign factor ±1 depending on that number. This approximation, along with the anti-commutation relations 
between the different 7 operators, allow us to move the operator 7^ from the far right of the above expression to 
just right of the operator j m> with only a minus sign appearing whenever m < i. The combination 7 m 7^ applied to 
(0| does not affect the state, and these two operators can be ignored. We therefore conclude that the effect of the 
occupation of a distant single-particle state is exponentially small and can be ignored. We now have to decide where 
to draw the line between states that we will consider and states that we will ignore. Since we are interested in lowest- 
order corrections, we shall consider only the effect of nearest-neighbour states, which will produce the lowest-order 
corrections. 

As a first step, let us consider the Rabi frequency of qubit i when driving only that qubit. Since under our 
approximation we only consider the effect of the occupation of the two neighboring sites, we only need to consider 
a three-qubit system (the algebra can be easily generalized to longer chains). It is straightforward to show that for 
driving qubit 1 we have 



I 71 (ci + cj) |0) = (0| 7271 (Si + c\) 7! |0) = (0| 7371 (&i + 4) 7^ |0) = («l 727371 (ci + 4) 7^ |0> 



^i(l) 



and a similar relation applies to driving qubit 3. 

(OI-72 (l-2c\ci 



C2 



For driving qubit 2 we have 
- t M«») = i -(2. 



(0| 7172 (l - 2c\c x ) (c 2 + 4) 7} |0) « ^ 2 (2) 

tI|o) 



1 



2:S 



(A 1 -A 2 )(A 2 -A 3 )(A 1 -A 3 ) 2 



(0| 7372 1 



2c|ci 



I 717372 (l ~ 2^ci) (c 2 + 4) 7 | 7 J |0) 



^2(2). 
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(A X -A 2 )(A 2 - A 3 )(Ai 



A3) 2 



(49) 

(50a) 
(50b) 

(50c) 
(50d) 



Therefore the Rabi frequencies of qubits 1 and 3 (which are at the end of the chain) do not depend on the state of 
the other qubits (this result holds even for longer chains and can be shown rigorously using Wick's theorem), while 
the Rabi frequency of qubit 2 shows some dependence on the states of the surrounding qubits, with a spread that is 
of the order of (J/SA) 4 . These results agree with the numerically derived results of reference 

We now turn to the Rabi frequencies when driving two qubits simultaneously in order to implement the SD-CNOT 
gate. We now have driving terms applied to two qubits. As a result, in addition to the expressions calculated above, 
we need to evaluate expressions of the form 



! 7m In ' ' ' 7p 



y<jj(l-2cjc,) (c 1+1 +4 +1 



3<i 



(51) 



As above, we ignore all the distant qubits. We consider a four-qubit chain with the driving applied to qubit 3 such 
that qubit 2 changes its state [i.e. i = 2 in equation (51)]. Here we take a chain of length four, to also incorporate 



the influence of the nearest-neighbour of the driven qubit. The relevant matrix elements are given by 
(0| 72 1] f 1 - 2 4%) + 4) 10} = lfe(3) 



(52a) 



3<2 



! 7i72 I] i 1 - 2d h) ( 5 3 + 4) 7i f |0) = ife(3) + 2^i(4) x 
i<2 



J12 Jf 3 J34 



(Ai - A 3 )(A 2 - A 3 )(A 2 - A 4 )(Ai - A 4 ) 
! 7472 1] i 1 - 2 4o) (c 3 + 4) t! |0) w ^ 2 (3) + 2^4(2)^4(3)^2(2) 



j<2 



" (0| 717472 J] ( X 2 4 a ( a 3 + 4) 7l7i f |0) w ^2(3) + 2^4(2)^4(3)^2(2), 

J<2 



(52b) 
(52c) 

(52d) 
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and similar expressions for the case where the operator 73 appears, e.g. 

- (0| 71737472 I] i 1 - 2d h) ( £ 3 + 4) lllhl 1°) = MQ> (53) 



J<2 

We therefore arrive at the main result of this section 



I 72 

j<2 



J<2 

I 7472 

J<2 



J] (l " 2ctq) (c 3 + 4) |0) = ^ 2 (3) (54a) 

J<2 

I 7172 J] ( X - 2 4^) ( £ 3 + 4) 7i |0> = lb, (3) ^1 + j (54b) 

J] (l 2c&) (c 3 + 4) 7! |0) = ^ 2 (3) ^l + O ^ (54c) 

- <0| 717472 II ( X - 2 %) ( 6 3 + 4) M |0) = ^(3) ^1 + (J^j j (54d) 

and similar expressions for the case where the operator 73 appears. 

The relative error in the CNOT gate will therefore be of order (J/SA) 2 . This result is in agreement with the 
numerical calculations in reference [55] where also a spread of (J/ 5 A) 2 is found in the matrix elements for flipping 
qubit i when driving qubit i + 1, depending on the states of the surrounding qubits. 

The scaling of the error with (J/SA) 2 is comparable to other schemes for implementing entangling two-qubit gates. 
The SD-CNOT gate has however the added advantage that the qubit level splittings do not need to be shifted. The 
shifting of energy levels gets increasingly complicated in multi-qubit systems, where the multitude of levels makes it 
harder to avoid crossing resonances that lead to unwanted evolution and entanglement. Also the freedom that the 
qubits do not need to be nearest-neighbour in frequency can make the design of a multi-qubit system simpler. 



VIII. CONCLUSION AND DISCUSSION 



We have theoretically analysed all the properties of the SD-CNOT gate that are important for its implementation 
in experiments. In the weak-driving limit, and in the absence of decoherence, the fidelity of the gate is 100%. For 
increasing driving strength the fidelity of the gate decreases due to the off-resonant driving of other transitions in 
the system. The off-resonant driving can induce spurious excitation of these transitions and AC-stark shifts of the 
levels, inducing an undesired single-qubit and (entangling) two-qubit evolution. We analysed the strength of these 
effects analytically, and in addition performed a full numerical simulation of the time evolution of the system. For the 
SD-CNOT gate without extra optimizations, the numerical simulations show a fidelity of 99% and 99.99% for gate 
speeds corresponding to 0.48J and 0.07J respectively. At high gate speeds the fidelity can be boosted by employing 
single-qubit correction pulses before and after the two-qubit driving. 

The gate performance was also analysed for multi-level qubits with weak anharmonicity of the energy splittings. 
As expected, it was found that the gate speed for a specified fidelity goes down for decreasing anharmonicity of the 
qubits. The energy shifts of the computational states due to the higher levels and leakage to these higher levels can 
both play an important role. We determined fidelities of 99% and 99.99% for gate speeds corresponding to 0.08d> and 
0.02(5 (for 5 sa 200 MHz). The gate speed can be increased by designing a relatively large coupling strength or low 
detuning between the qubits, or both. Another possible route towards increasing the gate speed and/or fidelity for 
weakly anharmonic qubits is using pulse shaping of the driving field. 

The scalability of the SD-CNOT gate was also studied. For a 1-dimensional array of qubits the error in the gate 
scales as (J/ 5 A) 2 , which is comparable to other schemes for implementing entangling two-qubit gates. The SD-CNOT 
gate has however the favorable properties that the qubit level splittings do not need to be shifted, and neither is it 
required that the qubits are nearest-neighbour in frequency. These properties can make a significant difference in 
the design of multi-qubit systems. We conclude that in terms of fidelity, gate speed and potential for scalability, the 
SD-CNOT gate provides a valuable addition to the existing variety of entangling gates. 
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Appendix A: Transformation to rotating frame 



We start from the Hamiltonian and driving term as given by equations ([I) and For this calculation we use 
the energy eigenbasis in the absence of the driving term [i.e. the states |0) , 1) , |2) and |3) in equations j9j]. The 
Hamiltonian without the driving term takes the form: 



ffn 



/- V /(A 1 +A 2 )2+4J2 



V o 



4 J 2 



V(Ai-A 2 )2 




4 J 2 



\ 



(Al) 



4J 2 / 



and the matrix elements of the driving term ffdrivc [equation (|3j)] are given by 



(1| F drl ve |0) 

(2|i? drivo |0) 

(3 1 -ffdrivc 1 2) 
(3 1 ffdrivo |1) 



(i\H, 



drive 



— a\ sin(9 + cos(wi + ipi) + a 2 cos(9_ cos(W + 992) 
+a% cos(? + cos(a;i + <pi) + a 2 sinf9_ cos(W + < f o 2 ) 
+ai sin(9 + cos(wt + + a 2 cos(9_ cos(W + <y2 2 ) 
+ai cos0+ cos(wi + </?i) — a 2 sin#_ cos(u>t + <^? 2 ) 

01 -ffdrivc \i) ■ 



(A2) 



All other matrix elements vanish. 

We now make a transformation to a convenient rotating frame. This means that instead of solving the Schrodinger 
equation for we solve an equation for 



w) = U{t) \ip) 



where 



U(t) 



iut/h 



/ -hu 





V huj 



(A3) 



(A4) 



The lab-frame state \ip) satisfies the Schrodinger equation: ih-^ = H In the rotating frame, the state 

should also satisfy the Schrodinger equation: ifi4z 
derived: 



Note that u is not the only possible choice that could serve our purpose here. Other transformations will lead to the 
same result. 

iji J. From these equations the transformation is easily 



1>) = H 
H = UHU^ - ihU 



dt 



With the total Hamiltonian given by ff 

H = UHqU^ - u 
= H -u 

( % /(A 1 +A 2 )^+4J2 
2 







= ffn 



= UHU^ - u. 
ffdrive, this results in 



(A5) 



V 

fft. 



\/ (Ax — A 2 ) 2 +4J 2 
2 










■ v /(A 1 -A 2 ) 2 +4./ 2 
2 








^(Ar+A 2 ) 2 +4J 2 
2 



(A6a) 



(A6b) 
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The matrix elements of i/drive can be conveniently split into H^n 
and W the time-dependent terms: 



V + W, where V contains the time-independent 



(1\V 


0} 


= e iut (l\V\0) 












— ( — n -, Qin n 

1 (_( I Oil! \J 

2 v 


+ e 


_|_ n o rns £ 

| Li ^ LUD 1/ 




(2\V 


o) 


i 

= - (+ai cos 




+ a.2 sin 6 




(3|y 


2} 


= — (+ai sin# 




+ a-2 cos 6 




(3|F 


1} 


= i (+ai cos 6 




— a2 sin 6 




(i\V 


li) 










h 1 w 


n\ 

/ 


— — — f {2,1 Sill $ 

2 V 


+e -t- 






(2| W 


0) 


= - (+ai cos 6 


4 


- a 2 sin 




(3| W 


2) 


= — (+ai sin 


+ 


a2 cos 9- 




(3| W 


1) 


= - (+ai cos 6 


+e' ;vi - 


d2 sin 




(i\W 


\j) 


= 0W)*- 









(A7) 



Transforming back to the lab-frame, the matrix elements in the energy eigenbasis of V and W are given by 



(1\V 


|0> = 


i 

2 


(— ai sin 9 + e~ 


-iipx 


+ 0,2 COS t 


?_e" 


-i V2 s 


e~ 




(2\V 


|0) = 


1 

2 


(+ai cos 6*+e~ 


-iipj. 


+ <Z2 sin ( 




-i V2 s 


e~ 




(3\V 


|2) = 


1 

2 


(+ai sin + e~ 


-iipx 


+ a 2 cos ( 


?_e" 


-i V2 s 


e~ 


iujt 


(3\V 


|1> = 


1 

2 


(+ai cos #+e~ 


-itpi 


— a2 sin ( 


?_e" 




e~ 


iujt 


{i\V 


li) = 


(j|V-|i>* 














(i\W 


li) = 


{' 


V\j)*. 















(A8) 

From these equations we see that V describes the co-rotating and W the counter-rotating field with respect to the 
Larmor precession of the system. 

The representation of V in the uncoupled-qubit basis can be straightforwardly calculated from the result of equa- 
tion (A8), however the result is a rather involved and large matrix. The main results of this appendix, i.e. the 



matrix elements in equation ( A7 ) , can be reproduced exactly by use of the terms H^ rivc and H^ iivc as described in 
subsection II A This approach is valid for the purpose of calculating transition strengths, but for other purposes it 
can be necessary to use the full representation V and W . 

Note that in the derivation in this appendix no approximations have been made: we have only performed transfor- 
mations and rearranged equations. 



Appendix B: Quantum mechanical description of the driving field 



In the other sections we have treated the driving field as a classical field. In this section we provide a description of 
the driven system where the driving field is treated quantum mechanically. With this description we verify some of 
the results on level shifts and gate errors obtained using the classical field description. Additionally, the description 
provided here can be used as a starting point to further investigate the strong-driving regime, maximum gate speed 
and the related gate errors. 

For this calculation we use the dressed-state picture. We consider a model composed of two qubits and a harmonic 
oscillator containing a large number of photons. In this picture the Hamiltonian becomes time independent and the 
frequency of the driving field now appears as hcub^b, where b and b' are the annihilation and creation operators of the 
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quantized field. The driving term dj cos(uit) in the original Hamiltonian (i.e. the Hamiltonian in which the driving 
field is treated as classical) is now replaced by the operator g^b + b>), where g t represents the coupling strength 
between the resonator and qubit i. 

We take the case where qubit 2 is the target qubit and where the 1 00) «-» 1 01} transition is to be darkened. In this 
case tpi — if2- Again, for simplicity in the expressions, we take fx=f2 = 0. Unlike in the weak-driving limit, we do 
not impose the relations Hlu = E\ — Eq and 02/01 = sin9 + / cos#_. Instead, we leave the coupling strengths between 
the oscillator and the two qubits as tunable parameters to be decided later. 

Assuming that the qubit frequencies Ai and A2 are much larger than the coupling and driving frequency scales, 
the energy-level structure breaks up into blocks of four energy levels in each block, with negligible effects caused by 
the coupling between the different blocks. We therefore consider the spaces spanned by the states |0, N + 1), |l,iV), 
|2, N) and |3, N — 1). In this space the Hamiltonian reads: 



H 



/ 

—C\ sin 9 + + C2 cos 9. 
Ci cos 9 + + C2 sin 9- 

V 



—C\ sin 9 + -t- 

(E 1 - Eq) - 


C\ cos 9 + — C2 sin 9. 



C2 cos 9- 

,q) — flLO 



Ci cos 9 + + C 2 sin 9_ 


(E 2 - Eq) - Huj 
C\ sin 9 + + C2 cos 6L 





Ci cos 9 + — C2 sin 6L 
C\ sin 9 + + C2 cos 6L 
{E x +E 2 - 2E ) - 2Huj 



(Bl) 



where Ci — y/Ngt, with N the number of photons in the resonator, and we have made the assumption that N is 
large, such that N w N + 1. 

One can now argue that the condition for a darkened transition is that the two corresponding energy levels must be 
degenerate: what corresponds to driven (Rabi) evolution in the previous picture can be regarded as Larmor precession 
in the current picture, and for Larmor precession to be suppressed the corresponding states have to be degenerate. 

If we make the approximation that the driving amplitudes are small compared to the detuning, the above Hamil- 
tonian can be approximated as 



H 



- Ci sin 6 + + C 2 cos 0_ 
-Ci sin 0+ + C 2 cos 9- \e x - E ) - huj 








(E 2 - E ) - hu 
C\ sin 9 + + C2 cos 6 






Ci sin 6+ + C 2 cos 6L 
{Ex +E 2 - 2E ) - 2hu 



(B2) 



and then it is obvious that the lowest two energy levels are degenerate for the choice of parameters (E\ — Eq) — Kio = 
and —C\ sin. + C2cos0_ = 0, which is just a different way of writing the condition derived in the semiclassical 



derivation [equation ( 13 1 



When we do not make the approximation of small driving amplitudes, using the full Hamiltonian in equation (Bl I, 
the matrix elements -ffo2> ^20 1 ^13 an d will mix all the quantum states to form the new, corrected energy 
eigenstates. If we think in terms of an order-by-order expansion of the effects of these matrix elements, we see that 
the matrix elements Hq2 and H2Q will cause mixing and repulsion between the states |0, N + 1) and |2, TV), while the 
matrix elements H13 and H31 will cause mixing and repulsion between the states |1, N) and |3, N — 1). To first order, 
the corrected quantum states are: 

|0,JV + 1)'= |0, N + 1) — Cl ™ 6 -EQ°-Z 6 ~ ' 2 ' < B3a ) 
|1, N)' = |1, N) - ClCO ;^ E ^Z e - 13, - 1) (B3b) 

|2, NY = |2, N) + C ^l C 2Z 9 ~ 1°' N + ( B3C ) 
|3, N - l)> = |3, N - 1) + II, N) , (B3d) 

and, to lowest order, the corrections to the energies (SEi = E[ — Ei) are given by 

SEq = -5E 2 = - V 1 cas0 + + C^0-)> 

E% — Eq — hU! 

8Ek = ~5E 3 = - (gic«g + -g2ring-) 2 (B4b) 

£/9 — En — Tun 
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With the above corrections, we have 



E[ 


-E' 


= E l 


-E + 


E' 3 


-E> 2 


= E 3 


-E 2 - 


E' 2 


-E' 


= E 2 


-E + 


E' 3 


-E[ 


= E 3 


-E x + 



4CiC 2 cos 9 + sin 6>_ 

E 2 — Eq — Tiui 
4CiC 2 cos + sin 0_ 

E 2 — Eq — huj 
2(CiC0S 0+ + C* 2 sin6L) 5 

E 2 — Eq — Tiuj 
2(CiCos6»+ - C 2 sin6L) 2 
E 2 — Eq — Too 



and 



(l,N\' H |0,JV + 1)' = -Ci sin (9+ + C 2 cos6L + (C x sin 6+ + C 2 cos 6L ) 



(CiCOB0+) 2 - (C 2 sin6>_) 2 

(E 2 -Eq- Twj) 2 



(3, N — 1|' H |2, N)' = d sin 6+ + C 2 cos 6L 



r „ ■ n n ACl COS 9+f - (C 2 Sm8_] 

+ (- Cl sm 0+ + C 2 cos 0_ ) (g2 + _ j go l M2 



(B5a) 
(B5b) 

(B5c) 

(B5d) 

(B6a) 
(B6b) 



The above equations can now be used to determine the modified conditions for performing the desired operation. 
One requirement for suppressing the |0, N + 1) «-» 1 1, iV) oscilla tions is that the parameters Ci, C 2 and <J are chosen 
such that the matrix element (1, N\' Hq |0, N + 1)' [equation (B6a)] vanishes. And in order to obtain full oscilla- 
tions between the states \2,N) and |3, TV — 1), the driving frequency u> must match the energy difference E' 3 — E' 2 
[equation (B5b)]. These equations can be solved iteratively starting from the weak-driving relations and calculating 
corrections. Note however that at some level of accuracy one would need to take into account the coupling between 
the different blocks in the dressed-state ladder. 

It is however clear from the above results that strong driving does not simply result in harmless modifications to the 
required driving frequency and the ratio between the two driving amplitudes. We identify three processes influencing 
the performance of t he ga te. T he firs t effect is related to the energy shifts in the transitions with splittings E' 2 ~E'q and 
E' 3 — E[ [equations (B5c) and (B5d)]. They contain a common shift of 2[(C± cos 9 + ) 2 + [C 2 sin 6L) 2 ]/(i? 2 — Eq — Twj). 
This shift represents the ac-Stark shift of the control qubit. For small 9 + and 6L we get 2C\/(E 2 — Eq — hoj), which is 
the analogue of the ac-Stark shift calculated in subsection IV B (an estimate for this shift using typical experimental 
parameters is provided there). 

Both the control qubit and target qubit transitions also get shifted with an amount that is opposite for the two 
associated transitions: ±4CiC 2 cos#+ sin#_/(i? 2 — Eq — Tuxi). What this result means is that during the operation of 

the CNOT gate, an effective term appears, which induces a (partial) CPHASE gate, and prevents us from 

obtaining an ideal CNOT gate. This error is the equivalent of the CPHASE-type error discussed in subsection |IV B| 

If we take J< Ai- A 2 and the case of strong driving (i.e. C\ ~ Ai — A 2 ) we find that the undesirable ai^a^ term 
has a coefficient that is of the order of J 2 /(Ai — A 2 ), which causes small errors on the timescale of a CNOT gate, 
which is given by J^ 1 . As discussed earlier, CPHASE-type errors can also be corrected with single-qubit gates. 

Lastly, when the matrix elements Hq 2 , H 2 q, H\ 3 and H 3 \ can not be neglected, it is not possible for any pair 
of energy levels to be degenerate, which is a general rule of quantum mechanics. As a result no transition can be 
completely darkened in this case, regardless of the choice of the parameters Ci, C 2 and u>. To first approximation this 
source of errors can be regarded as performing an additional single-qubit rotation, which is relatively easy to correct 
for in experiments. 

The results provided in this appendix can be used to perform a detailed numerical study for optimizing the different 
parameters in order to maximize the fidelity of the CNOT gate in the case of strong driving. However, we shall not 
do that here. 
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